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Abstract 



Spatiotemporal features and physics of vehicular traffic congestion occurring due 
to heavy freeway bottlenecks caused by bad weather conditions or accidents are found 
based on simulations in the framework of three-phase traffic theory. A model of a 
heavy bottleneck is presented. Under a continuous non-limited increase in bottleneck 
strength, i.e., when the average flow rate within a congested pattern allowed by the 
heavy bottleneck decreases continuously up to zero, the evolution of the traffic phases 
in congested traffic, synchronized flow and wide moving jams, is studied. It is found 
that at small enough flow rate within the congested pattern, the pattern exhibits a non- 
regular structure: a pinch region of synchronized flow within the pattern disappears 
and appears randomly over time; some of the wide moving jams merge randomly into 
a mega-wide moving jam (mega-jam). We show that the smaller the average flow 
rate allowed by the heavy bottleneck within the congestion, the less the mean pinch 
region width and the greater the mean mega-jam width; when the bottlenecks strength 
increases further, only the mega-jam survives and synchronized flow remains only 
within its downstream front separating free flow and congested traffic. Theoretical 
results presented explain empirical complexity of traffic congestion caused by bad 
weather conditions or accidents. 



1 Introduction 

The physics of freeway traffic congestion is one of the most quickly developed fields of 
complex spatiotemporal systems. In empirical observations, traffic breakdown (onset of 
congestion) in free flow occurs mostly at bottlenecks associated with, e.g., on- and off- 
ramps. In congested traffic, moving jams are observed [1-12]. A moving jam is a localized 
structure of great vehicle density, spatially limited by two jam fronts; the jam propagates 
upstream; within the jam vehicle speed is very low. 

Moving jams, which exhibit a characteristic jam feature [J] to propagate through bot- 
tlenecks while maintaining the mean velocity of the downstream jam front, are called wide 
moving jams [12]. A wide moving jam consists of alternations of regions in which traffic 
flow is interrupted and moving blanks [13, 14]; a moving blank is a blank between vehicles, 
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which moves upstream due to vehicle motion within the jam. The flow interruption can be 
used as a microscopic criterion for a wide moving jam associated with the jam definition 
[J] [13, 14]. 

In observations [12], traffic breakdown is associated with a local first-order phase tran- 
sition from free flow to synchronized flow (F— >S transition) at the bottleneck; synchronized 
flow [S] is defined as congested traffic that does not exhibit the feature [J]; in particular, 
the downstream front of synchronized flow is often fixed at the bottleneck. Wide moving 
jams can emerge spontaneously in synchronized flow (S^J transition) only, i.e., due to a 
sequence of F— >S— >J transitions. 

Moving jam emergence in synchronized flow leading to S— >J transitions is called the 
pinch effect occurring within an associated pinch region of synchronized flow (see Sect. 
12.2 of [12]): (i) The density increases and average speed decreases; narrow moving jams, 
which do not exhibit the feature [J], appear in the pinch region, (ii) These jams propagate 
upstream growing in their amplitude; as a result, S— J transitions occur, i.e., wide moving 
jam emerge. The upstream boundary of the pinch region is a road location at which a narrow 
moving jam has just transformed into a wide moving one. (iii) These locations can vary for 
different wide moving jams, i.e., the pinch region width depends on time. A congested 
traffic pattern that exhibits this regular structure is called a general pattern (GP). 

Earlier traffic flow theories and models reviewed in [1-10] cannot explain F— >S^J tran- 
sitions and the pinch effect (see a criticism of the theories in Ref. [12, 15]). For this reason, 
the author introduced a three-phase traffic theory (references in [12]) in which there are (i) 
the free flow, (ii) synchronized flow, and (iii) wide moving jam phases. The synchronized 
flow and wide moving jam phases associated with congested traffic are defined via the em- 
pirical definitions [S] and [J], respectively. The first three-phase traffic flow models showing 
the pinch effect are stochastic microscopic models [16, 17]. Later, other three-phase traffic 
flow models were developed [15, 18-22]. 

In general, the average flow rate g( con s) [j] within a congested traffic pattern upstream 
of a bottleneck is the smaller, the greater the bottleneck influence of traffic (called the 
bottleneck strength). Empirical data show that the flow rate g( con s) within GPs occurring at 
usual bottlenecks like on- and off-ramp bottlenecks, which is equal to the average flow rate 
in the pinch region q(P mch ), is approximately within a range 

g (con g ) = g ( P inch) = noo _ 17Q0 vehides/h/lane. (1) 

Features of GPs and other congested patterns occurring at usual bottlenecks determined 
by road infrastructure (on- and off-ramps, road gradients, etc.), for which condition £0) is 
valid, have already been studied in detail [12]. In contrast with the usual bottlenecks, due 
to bad weather conditions or accidents heavy bottlenecks can occur, which exhibit a much 
greater influence on traffic (greater bottleneck strength) that limits g( con s) to very small 
values, sometimes as low as zero. Features of traffic congestion at the heavy bottlenecks 
are unknown. 



1 The time interval of flow rate averaging is suggested to be considerably longer than time distances between 
any moving jams within a congested pattern. 
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As follows from recent empirical studies of sequences of wide moving jams [13, 14], 
the flow rate gr( blanks ) of low speed states associated with moving blanks within the jams is 
approximately within a range 

^(blanks) = 300 _ 60Q veh i c i es / h /i ane . (2 ) 

The flow rate between the jams associated with non-interrupted flows in the jams' outflows 
q^J t is greater than q(P mch ) (TJJ), i.e., qj^ t is considerably greater than gr( blanks ) Now we 
assume that due to bad weather conditions or an accident a heavy bottleneck occurs with a 
great strength for which 

g(cong) ^(blanks) _ ^) 

In this case, q^ t must reduce to g( blanks ) (ffj, i.e., the difference between flows within and 
outside wide moving jams disappears. As a result, all wide moving jams should merge 
into one mega-wide moving jam (mega jam for short). Thus already from this qualitative 
consideration, we can expect interesting physical phenomena associated with complexity 
of traffic congestion at heavy bottlenecks. 

In this letter, we reveal these features and compare them with measured congested pat- 
terns at heavy bottlenecks caused by bad weather conditions or accidents. 



2 A Theory of Traffic Congestion at Heavy Bottlenecks 

For a theoretical analysis of traffic congestion at heavy bottlenecks, we use a stochastic 
three-phase traffic flow model of a two-lane road [16,23] (see Appendix) with the following 
heavy bottleneck model. 

We suggest that there is a section of the road within which due to an accident or bad 
weather conditions drivers should increase a safety time gap r( safc ) to the preceding vehicle. 
We simulate this effect by an increase in r( safc ) to some r( safc ) = Tb > 1 sec in comparison 
with r( safe ) = 1 sec used in the model for vehicles moving outside this section. In according 
with Eq. ([TBI of the model (Appendix), r( safc ) determines a safe speed, which should not 
be exceeded by a driver; otherwise, the driver decelerates. As a result, within this section 
drivers move with time headways that are approximately equal to r^ safe ^ = Tb. Therefore 
the section with longer r( safc ) = Tb acts as a bottleneck on the road. In the bottleneck 
model suggested here, each chosen value Tb defines a specific bottleneck: the strength of 
this bottleneck is the greater, the longer Tb; in turn, the longer Tb, the smaller g( con s), 
i.e., the greater the flow rate limitation within congestion caused by the bottleneck. This 
model feature allows us to simulate a heavy bottleneck caused by bad weather conditions 
or accidents leading to a long enough Tb within the bottleneck 0. 

2 Under bad weather conditions, a road section with a long value Tb can be caused, e.g., by a poor view due 
to fog on the section or a much longer deceleration way needed by snow and ice on the section. If an accident 
occurs on a road, a road section with a long value Tb can be caused, e.g., by much narrower lane widths allowed 
for driving on the road section; the same effect can occur under heavy roadworks. 



4 



Boris S. Kerner 



When Tb is chosen to be not very long, then at a great enough flow rate q\ n in free flow 
upstream of the bottleneck, firstly an F^S transition occurs spontaneously at the bottleneck. 
Then the pinch region with a relatively great flow rate g( con s) = q(P inch ) is formed. At 
the pinch region upstream boundary wide moving jams emerge. Thus we found known 
phenomena of regular GP formation (Figs. Q] (a), [2] and [3] (a, b)) [12]. The pinch region 
width L (pinch) (t) changes over time between about 1 and 2 km (Fig. H (a)). L (pinch) is 
defined as the distance between the upstream boundary of the bottleneck (x = 16 km) and 
the road location upstream at which a wide moving jam has just been identified through 
the use of the jam microscopic criterion. There is also a nearly constant frequency of 
£ (pinch) ^ oscillations associated with the maximum in the Fourier spectrum (Fig. H](b)). 
Speed autocorrelation functions and associated Fourier spectra of speed time-dependencies 
at shorter Tb show regular character of wide moving jam propagation (Fig.[3](c, d)). 

When Tb increases and therefore g( con s) = q(P mch ) decreases (Fig. [2]), and Tb remains 
a relatively small value (1.8 < Tb < 3 sec), then as for other bottleneck types, we found 
the following known GP features [12]: the smaller the flow rate g(P mch ), the greater the fre- 
quency of narrow moving jam emergence within the pinch region, the lower the maximum 
speed between wide moving jams upstream of the pinch region, the smaller the mean pinch 
region width L^n h) (Figs. ffl(b) andi(c, d, i)). This can also be seen from a comparison 
of time-dependences of average speed within the region of wide moving jams for different 
T B (Fig.[3](a,e)). 

Qualitatively other phenomena are found when Tb further increases (Tb > 3 sec) and 
the average flow rate g( con §) decreases considerably (Figs.[T](c, d) and[2]). Firstly, L(P mch ) (t) 
becomes a non-regular time-function (Fig.|4](e, g)) whose Fourier spectrum is broader, the 
longer Tb (Fig.|4](f, h)); there are random time intervals when the pinch region disappears, 
i.e., £(P mch ) = o (Fig. 01(e)). This means that there are time instants at which there is 
no pinch region and wide moving jams emerge directly at the upstream boundary of the 
bottleneck, whereas for other time intervals the pinch region appears again (Fig. [4] (e)). 
Consequently, L^^t^ decreases continuously, when T B increases 
reaches almost zero for Tb = 30 sec: for such a heavy bottleneck, the pinch region does 
not exist. 

When Tb increases and the pinch region disappears during some random time intervals, 
the GP behaviour qualitatively changes: the time-functions of average speed become non- 
regular ones (Fig. [3](i, m)). We found that under condition ©, the pinch region disappears 
and wide moving jams merge into a mega-jam. In particular, for Tb = 30 sec no separated 
wide moving jams can be distinguished: instead of the sequence of wide moving jams, a 
mega-jam occurs only. These effects can also be seen from speed autocorrelation functions 
and associated Fourier spectra of the average speed time-dependences (Fig. [3](k, 1, o, p)). 
The mega-jam is found at those long values of Tb at which the pinch region disappears. 
Simulations show that the mega-jam consists of alternations of flow interruption and low 
speed states associated with moving blanks. Thus there are no continuous flows within the 
congested pattern any more with the one exception of the downstream front of synchronized 
flow that separates free flow downstream and the mega-jam upstream of the front. Synchro- 
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Figure 1: Simulated speed in space and time in the left (left) and right (right) road lanes 
at different T B : T B = 1.8 (a), 2.4 (b), 12 (c), 30 sec (d). q in = 1946 vehicles/h/lane. The 
upstream boundary of bottleneck region of the length 300 m is at x = 16 km; the maximum 
speed wthin this region is 60 km/h. Resulting values of q( con & = 1546 (a), 1114 (b), 440 
(c), 217 vehicles/h/lane (d). 
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Figure 2: Simulated average flow rate q( con s) within congested patterns shown in Fig. [T]as 
a function of Tb . The averaging time interval for q( con & is 60 min. 



nized flow remains only within this front. Thus in accordance with all known empirical data 
and theoretical results presented, the condition 

g(cong) > ^(blanks) ^ 

is a necessary condition for GP existence at a bottleneck. 

For explanations of these results, we recall that an S— >J transition is a first-order phase 
transition [12], i.e., it is characterized by a random time delay Tgj. Tgj includes a random 
time delay of spontaneous nucleation of a narrow moving jam (nucleus for the S— J tran- 
sition) and its random growth within the pinch region until the jam transforms into a wide 
moving jam at the upstream boundary of the pinch region. Thus the smaller Tsj, the shorter 
.L (pinch). The random character of Tsj explains a time-dependence of l/ pmch )(i) (Fig.@](c, 
e)). We found that the greater Tb, the greater the density within the pinch region and the 
smaller the mean random time delay 7^™ ean ) Thus the greater Tb, the greater the proba- 
bility for the occurrence of a negligibly small Tgj and wide moving jam emergence directly 
upstream of the synchronized flow front; in the latter case 

L ( P inch) = o This explains also 

why Lmean^ decreases up to zero, when an increase in Tb causes a decrease in 7 , g (™ can ) U p 
to zero, i.e., when all wide moving jams emerge almost directly upstream of the bottleneck 
and the pinch region disappears. 

When g( con s) becomes zero, because behind a road location the road is closed, the 
mega-jam transforms into a queue of motionless vehicles, which therefore is not associated 
with vehicular traffic. Nevertheless, there is a link between the queue and the mega-jam. If 
at a time instant we allow several vehicles to escape from this queue, then simulations show 
that motion of these vehicles results in wide moving jam occurrence: the downstream front 
of the jam separates moving vehicles escaping from the queue and vehicles standing within 
the queue upstream of the front. When the number of vehicles escaping from the initial 
queue decreases, the downstream jam front transforms into a moving blank(s) subsequently 
covered by vehicles standing in the queue. When a vehicle per a long enough time interval 
is allowed to escape from the mega-jam, as simulations show, a sequence of such moving 
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Figure 3: Characteristics of congested patterns shown in Fig. Q] related to location 10 km. 
Time-functions of speed (a, e, i, m), speed correlations (c, g, k, o), associated Fourier 
spectra (d, h, 1, p), and flow rate (b, f, j, n) for different T B = 1.8 (a-d), 2.4 (e-h), 12 (i4), 
30 sec (m-p). 1-min average data. 
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Figure 4: Simulations of pinch effect: L( pinch ) (t) (a, c, e, g), associated Fourier spectra (b, 
d, f, h) for congested patterns shown in Fig. ED L^an h) (g (cong) ) (i). In (a-h) T B = 1.8 (a, 
b), 2.4 (c, d), 12 (e, f), 30 sec (g, h). 
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blanks within this jam occur; these moving blanks exhibit, however, dynamics specifically 
associated with the mega-jam. 

Thus based on a study of three-phase traffic flow model we found that the complexity 
of traffic congestion caused by bad weather conditions or accidents, when a very heavy 
bottleneck appears on a road, is associated with the phenomenon of random disappearance 
and appearance of the pinch region over time that is accompanied by the random merger of 
some wide moving jams into a mega-jam. When the bottleneck strength increases further 
strongly, the pinch region disappears and only the mega-jam survives and synchronized flow 
remains only within its downstream front separating free flow and congested traffic. These 
phenomena explained by random nucleation of moving jams in congested traffic reveal the 
evolution of the traffic phases when heavy bottlenecks occur in highway traffic. 

3 Discussion. Comparison with Empirical Results 

To compare the above theory of traffic congestion with empirical congested patterns, one 
should have measured data for traffic congestion at a bottleneck, whose strength should be 
manually continuously changeable from the one associated with usual bottlenecks like on- 
and off-ramps to great bottleneck strengths associated with heavy bottlenecks caused by 
bad weather conditions or accidents. Unfortunately, such measured data is not available. 
However, we can compare the theory with measured data related to two limiting cases: (i) 
Traffic congestion at an usual on-ramp bottleneck (Fig. [5]). (ii) Traffic congestion at a heavy 
bottleneck caused by bad weather conditions - snow and ice on a road (Fig.©. 

The speed distributions with congested patterns found in simulations (Fig. [T](a, b)) for 
the range of Tb = 1.6 — 2.4 sec associated with the flow rate range 

^(cong) = g ( P inch) = 112Q _ 18Q0 vehic i es /h/l a ne (5) 

are qualitatively the same as those in an empirical GP shown in Fig. [5] and in all other 
known empirical GPs [12]. Moreover, quantitative values of empirical flow rates q(P mch ) 
dD are approximately associated with the theoretical result ©. In particular, within the 
pinch region of the GP shown in Fig. [5] (a), the flow rate averaged between 7:00 and 8:00 
and across the road is g(P mch ) = 1200 vehicles/h/lane. 

In contrast, for Tb > 6 sec associated with the average flow rate 

g (con g ) < 625 vehicles/h/lane (6) 

the theoretical speed patterns (Fig. Q] (c, d)) are qualitatively the same as those we found 
in empirical data measured on many various days (and years) on different freeways, when 
very heavy bottlenecks caused by bad weather conditions or accidents occur for which 
empirical flow rates g( con s) < 700 vehicles/h/lane. In this case, rather than regular structure 
of traffic congestion of GPs (Figs. Q] (a, b) and [5]>, both empirical and theoretical traffic 
congested patterns (Figs. Q] (c, d) and [6]) exhibit non-regular spatiotemporal structure of 
congestion in which no sequences of wide moving jams can be distinguished. In addition, 
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Figure 5: Usual empirical traffic congestion - a general pattern (GP) at on-ramp bottleneck: 
(a) Speed in space and time for the GP. (b) Speed in different lanes within a sequence of 
wide moving jams associated with (a) at location 1.5 km. 1-min average data measured on 
the freeway A5-South in Germany dated April 15, 1996. Arangement of detectors is shown 
in Fig. 2.1 of the book [12]. 
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the mentioned empirical flow rates g( con §) < 700 vehicles/h/lane found within the empirical 
non-regular congested patterns are associated with the theoretical result (©. In turn, these 
theoretical and empirical results for q( con s) correspond to empirical flow rates found within 
wide moving jams associated with moving blanks (T2]>. 

An example of such an empirical congested traffic pattern is shown in Fig. [6] Indeed, 
rather than the regular structure of congestion within the GP (Fig. [5]), in measured data 
associated with bad weather conditions a non-regular spatiotemporal structure of congestion 
is observed (Fig. [6]). A heavy bottleneck appears on February 02, 2006 between locations 
4.07 and 3.02 km due to snow and ice. Upstream of the bottleneck, very low speed and flow 
rate patterns (x < 3.02 km in Fig.[6](b)) are observed. Downstream of the bottleneck (x = 
4.07 km) vehicles have escaped from the congestion (speed is high), however, the flow rate 
is very small because the bottleneck reduced the average flow rate within the congestion 
strongly. For example, at x = 3.02 km, the flow rate q( con & averaged between 7:00 and 
7:40 and across the road is only 513 vehicles/h/lane that corresponds to the theoretical result 
©. In contrast with the GP in Fig.|5](b) (as with other empirical GPs [12]), within traffic 
congestion shown in Fig. [6] (b) non-regular low speed patterns are observed in which a 
sequence of wide moving jams cannot be distinguished (x < 3.02 km). This conclusion 
is regardless of the flow rates to on- and off-ramps, which in the data set lead to a greater 
average flow rate downstream of the freeway intersection (Fig. [6] (a, b)). 

If speeds in different lanes are compared (Fig. [6] (c)), we find even more non-regular 
speed time-dependencies: whereas no vehicles pass a detector in one of the lanes, i.e., the 
speed is zero, during the same time interval the average speed in other lanes can be higher 
than zero. This explains why in time-dependences of speeds averaged across the road (Fig. [6] 
(b)), this very low speed is seldom equal to exactly zero. 

I thank Sergey Klenov for help in simulations, Ines Maiwald-Hiller for help in data 
analysis, Andreas Hiller, Gerhard Nocker, Hubert Rehborn and Olivia Brickley for sugges- 
tions. 

A Stochastic Microscopic Three-Phase Traffic Flow Model 

The stochastic microscopic three-phase traffic flow model for two-lane road used in all 
simulations presented above reads as follows [16,23]: 



v n+ i = max(0, min(ff ree , v n+1 + £ n , v n + ar, u S) n)), 



(7) 



X n +1 = X n + V n+ lT, 



(8) 



v n+ i = max(0,mm(u fre e,fc,n,fs,n)) 



(9) 
(10) 




v n + A n at g n < G 
v n + a n r at g n > G 



n 



(11) 
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Figure 6: Empirical structure of congestion caused by snow and ice: (a) Schema of road 
detector arrangement on a section of the freeway A3-North in Germany near the intersection 
"Hanau". (b) Average speed and flow rate across the freeway at different locations, (c) 
Speed in different lanes at two locations. 1-min average data. 
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where 

A n = max(-6 n r, min(a n r, v^ n - v n )), (12) 

where index n corresponds to the discrete time t = nr, n = 0, 1, 2, ...; r is time step; v n 
is the vehicle speed at time step the vehicle co-ordinate; v n is the speed calculated 

without a noise component £ n ; Vf iee is the maximum speed in free flow that is constant; a is 
the maximum acceleration; g n = x# n — x n — d is the space gap; the lower index i marks 
functions and values related to the preceding vehicle; v Sjn is a safe speed; all vehicles have 
the same length d; a n > 0, b n > (see below); G n is a synchronization gap: 

G n = G(v n ,v ttn ), (13) 

where the function G(u, w) is chosen as 

G(u, w) = max(0, km + (f>oa~ 1 u(u — w)), (14) 

k > 1 and 4>q are constants. 

The safe speed v s n in (O is chosen in the form 

v s , n = min (v n saie \g n /r + vf>), (15) 

where Vn ^ = v^ aie ^ (g n ,V£ jn ) [24] is a solution of the Gipps-equation [25] 

u{ safc) T ( S afe) + X<J = ^ + X d (v e , n ), (16) 

where Xd(ti) is the distance travelled by the vehicle with an initial speed u at a time- 
fa) 

independent deceleration b until it comes to a stop; v e is an "anticipation" speed of the 
preceding vehicle at the next time step (formula (16.48) of [12]). 

The noise component £ n in |7]) that simulates random deceleration and acceleration is 
applied depending on whether the vehicle decelerates or accelerates, or else maintains its 
speed: 

r ifs n+1 = -i 

Hn = { U ifS n+1 = l (17) 
[ if 5 n+1 = 0, 

where S in (fTTT ) denotes the state of motion (<S n+ i = — 1 represents deceleration, S n+ i = 1 
acceleration, and 5 n+ i = motion at nearly constant speed) 

{-1 if v n+ i <v n -5 
1 if v n+ i >v n + 5 (18) 
otherwise, 



5 is a constant (5 <C ar), 



£ a = arG(p a - r), 
£b = aT@(p h - r), 



(19) 
(20) 
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where p a and p^, are probabilities of random acceleration and deceleration, respectively; 
r = rand(0, 1), &(z) = at z < and G(z) = 1 at z > 0. 

To simulate driver time delays in acceleration or deceleration, a n and b n are taken as 
stochastic functions 

a n = a&(P Q -r 1 ) (21) 

6 R = oe(Pi-n) J (22) 
Po = \i ifS n = l, (23) 

Px = ( Pl r , (24) 

{ P2(Vn) if 5 n = -1, 

where speed functions for probabilities Po(v n ) and P2(v n ) are considered in [12]; p\ is 
constant; r\ = rand(0, 1). 

The following incentive conditions for lane changing from the right lane to the left lane 
(R — > L) and a return change from the left lane to the right lane (L — > i?) have been used 
in the model: 

R^L: v + > t^ n + 5i and u n > (25) 

L — > i? : u+ > + 5i or t;+ > <y n + <5i. (26) 
The security conditions for lane changing are given by the inequalities: 

5+ > min(t; n T, (27) 

5~ > min(?;~T, G~), (28) 

where 

G+ = G(« nj v+), (29) 

G- = G(v-,v n ), (30) 

the function G(u, w) is given by (fT4l : the speed v£ or the speed ve >n in (|27T ). (T28T ) is set to 
oo if the space gap g+ or the space gap g n exceeds a given look-ahead distance L a , respec- 
tively; superscripts + and — in variables, parameters, and functions denote the preceding 
vehicle and the trailing vehicle in the "target" (neighboring) lane, respectively The target 
lane is the lane into which the vehicle wants to change. If the conditions (|25T>-(|27T) are 
satisfied, the vehicle changes lanes with probability p c . p c , 5\ (8\ > 0), L a are constants. 
Explanations of the physics of this traffic flow model are given in the book [12]. 
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